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Abstract 

We apply recent results in the theory of PDE, specifically in problems 
with two different time scales, on Einstein's equations near their Newtonian 
limit. The results imply a justification to Postnewtonian approximations when 
initialization procedures to different orders are made on the initial data. We 
determine up to what order initialization is needed in order to detect the 
contribution to the quadrupole moment due to the slow motion of a massive 
body as distinct from initial data contributions to fast solutions and prove 
that such initialization is compatible with the constraint equations. Using 
the results mentioned the first Postnewtonian equations and their solutions 
in terms of Green functions are presented in order to indicate how to proceed 
in calculations with this approach. 



1 Introduction 



In recent papers there has been given a rigorous justification of the Newtonian 
limit approximation in General Relativity. In one case, assuming symmetric 
hyperbolic equations for the matter (including appropriate boundary conditions for 
it), and in the other, assuming Vlasov type matter, there has been shown that 
given a Newtonian solution, there exists a nearby general relativistic solution for 
a time intervall which is independent on the limiting parameter. The proofs of 
these results relay on ideas pioneered by H-0 Kreiss on dynamical systems with 
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different time scales which in turn relay on energy estimates for symmetric hyperbolic 
systems. 

The above mentioned results are based on a initialization procedure, that is, the 
proximity of full relativistic solutions to the Newtonian ones is obtained by choosing 
the initial data in a very special way, basically ensuring that the time derivatives 
of the solutions at the initial surface stay bounded on that limit. This initialization 
procedure is only needed to a finite order on the limiting parameter (which is taken 
to be one over the speed of light) implying that fast behaving parts of solutions stay 
under control order by order along evolution. 

In this paper we advance further into this problem by showing two things: 

First that the General Relativity equations can be cast into a form on which 
the standard theory of different time scales applies (c.f. [§], [Q, and P). This 
substantially improves the results in [|I| , for there there were singular terms outside 
the principal part which had to be dealt with in a very involved way. This gives 
further information on the behavior of the fast part of the solutions, that is the part 
that comes from the failure to initialize data to all orders. In particular, it gives the 
equations that the leading order fast behaving parts of the solutions satisfies and 
so its bulk behavior. This result should be important in studying the scattering of 
gravitational waves on slow varying sources. 

To be able to apply these standard results (estimates) we shall assume that 
sources for Einstein's equations also satisfy symmetric hyperbolic equations which 
are regular in the limit and, if boundary conditions are needed to deal with them, 
that they are of such nature as to allow for the estimates to hold. Admittedly there 
are some problems to grant this to hold for some specific cases, but that is a problem 
on our understanding of the description of normal matter, and are not very much 
related to the dynamics of the gravitional degrees of freedom, so we do not address 
this issue further. 

Second that the initialization procedure can be done to higher enough orders so 
that the quadrupole formula should follow. The initialization procedure implies some 
relations between initial data besides the one implied by the constraint equations, 
and so the above mentioned results implies that showing that certain elliptic systems 
of equations have solutions with the appropriate asymptotic behavior. The result is 
only an argument, for the method used does not allow us to have estimates, and so 
to control evolution, all the way up to some portion of future null infinity, the region 
where the quadrupole formula should hold. A rigorous result for that issue should 
follow by studying evolution either along null cones, as considered by Winicour, , 
or along asymptotically null surfaces as studied by Friedrich, 0. 

The plan of the paper is as follows: 

In the second Section we introduce the general theory of different time scales 
systems and quote the relevant Theorem. We then follow |I|] and cast Einstein's 
equations as a symmetric-hyperbolic-elliptic system. By assigning units to the fields 
appearing in those equations one can determine how the parameter e = -, where c 
is the velocity of light, appears in the system and in this way the solutions become 
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a one parameter family of solutions to Einstein's equations. A further rescaling of 
the fields is needed so that the structure of the equations is in the form assumed in 
the hypotesis of the mentioned Theorem. 

In Section 3 we calculate the fiux of energy coming from the slow part of the 
solution (i.e the contribution to the fiux of energy due to the slow motion of a 
massive body) and the flux of energy coming from the initial data without sources, 
that is the contribution from fast solutions. We conclude that one should initiahze 
up to order three in e in order to isolate the contribution to the energy fiux coming 
from the slow part of the solution. Wc further prove a Lemma stating that such 
initialization is consistent with the constraint equations and so that solutions with 
these characteristics exists. 

In an Appendix we look at the solution to the initialized data to order e and 
obtain an explicit solution in terms of the shift vector N"- and the source S"^. This is 
done in order to illustrate how the general setting of the theory of different time scales 
produces the correct Postncwtonian equations to that order and how to proceed if 
one whishes to compute higher order corrections. 



2 Preliminaries 

The study of systems with different time scales reduces to the study of systems of 
partial differential equations which are singular in the limit of one of this time scales 
going to zero. The study of such a limit distinguishes between different classes of 
solutions according to their limiting behavior. Solutions which behave smoothly on 
the limit are called slow, they move according to the time scale which remains finite, 
those which do not have a well defined limit are called fast, they have dependence 
on the time scale we are setting to zero. As an example we consider the following 
system: 

(2.1) Ut 

(2.2) vt 

Since the system is linear we can consider its Fourier modes, with the ansatz 
we obtain 

It is clear that the first term in this expression is a fast solution, while the 

second term is a slow one. Notice that for a solution to be slow it is necessary 
that its time derivative stay bounded, in particular at t = 0. One of the main 
results of the theory of different time scales is that, under some circunstances, this 
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condition is also sufficient. In particular this allow us to pick slow solutions from 
initial data requiring that a number of consecutive time derivatives at the initial 
time are bounded, this procedure is called initialization. Another interesting output 
of the theory is that it gives precise information on the evolution of the fast and 
slow parts. 

We turn now to the general setting and quote the main results of the theory in 
the following Theorem. 

Theorem 2.1 ([|], see also Let 

(2,3) A%isu^)l-u' = (i/r, + Ai.,(As))a.^ + B.iu\s), 

be a symmetric-hyperbolic system containing the small parameter e; i.e., let the 
matrices K"^ , A\ and be symmetric and positive definite. Assume that the 
matrices A1 and the vector B are continuous in e uniformly in u for bounded u, 
and are in u uniformly in e with s > Sq = [n/2] + 2 and that the initial data 
u{0,x,e) = {qab{x,0),pab{x,0),r'^f,c{x,0)) lies in (RT") . Then the solution u(t,x,e) 
exists for a time T independent of e. Moreover if the initial data for \(2.3]\ is 

u{0, X, e) = uo{x) + e ui{x) + 0{e'^) 

where 

K^daUo = 0, 

then 

u{t, X, e) = x) + e[t;°(t, x) + u{t, x, e)] + C(e^), 
with and u satisfying 



A°(OX + A!J(M°,0)9aM° + K'^9at;° = Bo 

(2.4) K'^daU^ = 

u^{0,x) = Uo{x) 

and 

[A''{0)+eu''A'l{0)]ut + A1{u^,0)daU + ^K''daU = B, 

(2.5) u{0,x) = ui{x) -v^{0,x) 
respectively 



Bo = 5(t,x,M°,0) 

Bi = Be{t,x,u°,0) + {v'^ + u)Bu{t,x,u'^) 

- (t;° + m)A;J(m°, 0)9aM° - A%0)v^ - A"(m°, 0)dav'', 

and f ° chosen in such a way that K"-davP = . 
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The solutions of |(2.5)| are fast solutions depending on the choice of the initial 



data. There exists further local (in T") refinements of these results, that give us 
information about the behaviour of these fast solutions, see for example [|lOl, fll 



and 0. These results should be an useful tool for the understanding of problems 
like the scattering of gravitational waves by slow moving bodies. 

If we would now want that the fast solution appears in the second order of e, 
and would have to satisfy equations similar to K2.4)| and u one similar to |(2.5 



(with different sources), so that Ut and Uu remain bounded when e — ^ at t = 0. 
Summarizing we can suppress the fast part up to 0{e"'), demanding that n-time 
derivatives stay bounded when e — > O.This procedure is called initialization, for 
more details see and 0. Note that in the case of General Relativity the choice 
of initial data for the initialization procedure is not trivial for it must be consistent 
with the constraint equations. 

In order to give a formulation of the Postnewtonian limit in General Relativity as 
an initialization procedure, we use the variables given in such that the Einstein 
equations can be put as a symmetric hyperbolic-elliptic system. Defining the lapse 
function as 

(2.6) ^ 



with U being the Newtonian potential and 

(2.7) r'^'c = ^(^cg"' - lq^'qe,d,q'') ^ ^^^^iV^Q^') 

(2.8) p^' = 4vr"^ 

where dc is the derivative operator associated to a flat three-metric e"'' and qab is 
the three metric induced on the hypersurfaces t = const.. The evolution equations 
become 

^ab ^ -Se^q^'^NU + 2e\l-4e^U)N\^q°-^p-p''^) 

(2.9) -2e^q^ND^''N^^ + 2e^q''^D^N\ 



p 



-q-^N^^-iq'"^ddr''\ - 2q<''dcr^'^'^ d) - 2q^^q''\/\U - p) 

e 

ah I ^2 7\rcQ nb , ^ jpab / ^ de de 



(2.10) +25"" + e'N^dcp"" + £F^\e, r^^^p"^, dJJ, d^N' 



r"^ = -^\d,p'''-25'fddp'^'') 
qA e 

+ -{q'''d,ddN''-q''^^d,ddN''^} 
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qi e e 
-2eNq''''N'^dcddU 



(2.11) 



+eF 



ab 



ab ^ab 



CI 



where A^" is the shift, Da is the covariant derivative associated with qab, ^ = q"'''dadb, 



?ab 
?ab 



?abi 



e, r 



ab ^ab 



F c{£, r c,p , dcU, U, dcN"'), and all F's in the equations are smooth point- 
wise functions of all their arguments and the quadratic terms in r°'^c and p"* are the 
highest powers the factors will appear in and are 0{e'^). 
In these variables, the constraint equations become 



(2.12) 
(2.13) 



AU - -ed.r"'^ = P + e'F{e, r^\, dJJ) 



- 2dcp'"' = Ar + e'F\e, r'^",, dJJ) 



(2.14) 



„ab 



ab\ 



The system |(2.10)| - [[2.11)| is symmetric hyperbolic, and has an energy estimate 
finite for fixed e . In order to have an energy estimate bounded in the limit when 
e ^ the authors in [|l|, choose a gauge (i.e. a selection of lapse and shift) that 
makes r°-''b = and r"-^b = 0. As a matter of fact one could relax this gauge choice 



to r°'^b = C>{e^) and r°'''b = 0{e^) where /? is a positive and arbitrary number that 
guarantees that the energy estimates stay bounded even when e 0. In this paper 
we choose (3 = 3, the reason of this choice depends on the initialization procedure 
and it shall become clear in the next Section. 

This choice implies the following equations for A^": 



(2.15) 



4(7" + d'NU) + 2e\''^N''dcddU - e'G' 







(2.16) 



DdN'' 



-2NU, 



combining equations (2.15) and (2.16) we get 



(2.17) d.D^N'' + d^D^N" - + 2e^q''^N^d^ddU - e^G^ 

and the constraint equation |(2.12)| becomes 



(2.18) 



Af/ 



p + e^F{e,r'^\y\d,U). 



Thus the system given by [f2l0]l , fLVTl , fHM , fHM , K^Tt]] and [[2J8 



con- 



stitutes a symmetric hyperbolic-elliptic system, the hyperbolic and elliptic part are 
given by the equations |(2.10)| , |(2.11)| and |(2.17)| , |(2.18)| respectively. 
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By inspection one can see that the source Bi{u'', e) in |(2.10)| and |(2.11)| takes the 



M^, e), where Cij is a constant matrix. Since we want to 



form e) = + F, 

apply Theorem \2.1\ on this system, the source has to be smooth in e, thus we rescale 
the variables^ in the following way. 



-.21 ab 



(2.19) ff = erf p"^ = ep"" and q"" = e"^ + e^h" 

Since F""^ and F^^ have quadratic terms in the variables rf' and p"'^ , we obtain a 
system such as |(2.3)| with a source term depending smoothly on e, namely: 



•ab 



(2.20) 



-Se'^q'^^NU + 2e{l - Ae^U)N^q^{q''^p 



p 



ab\ 



= -qlN^}-[q'^dQ^-ab^ _ 2q<''d,f'''>'^d) - 2q^p''\AU - p) 
(2.21) +268"^ + e^N'dcp"^ + eF''\e, f'^\, f\ dJJ, d^N^) 



f'-'c = -^-{dcP^''-26i^d,p'^') + 
qi 6 

+q''^d,ddN^ - q'^^^dcddN''^ 
1 

-2e^Nq''^N'^d^ddU 
(2.22) +e^F-\{e, f"^,p'^^ d,U, U, a,iV"). 

In the same gauge as before, the constraint equations and the equation for the 
shift vector become 



(2.23) AU = p + e'F{e, f'K, d,U) 

(2.24) - 28^"'' = AeJ" + e^F^ie, f'^^,p''^ d^U) 



(2-25) r-^ = 7^^9,(7^5'^^), 



1 

2E\/q' 
and 

(2.26) dcD^N'' + d^D^N" - 4J^ + 3eF^ + e^{ANq^^N'^dcddU - F''\ 



similar rescaling is used in 
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From now on, to simplify the notation, we will work without bars on the variables 
and r°*c- As stated in the Introduction we assume that the above system is part of 
a bigger symmetric hyperbolic system which includes the equations for the sources, 
they are assumed to be regular in e and of such a type that global estimates can be 
obtained. 

The existence of smooth one parameter family of solutions is guaranteed because 
they prove that energy estimates corresponding to the hyperbolic part of the system 
stay bounded when e ^ via a particular gauge choice, the elliptic variables satisfy 
a Garding estimate in terms of the hyperbolic ones and the boundedness in time of 
the solutions are guaranteed via the initialization procedure. 

Another interesting problem which can be tackled with these techniques is whe- 
ther for each given fast Postnewtonian solution (i.e. a fixed slow background solution 
and a highly oscillating part), there exists a solution to the full Einstein equations in 
the same intervall of existence and that it remains close for that whole intervall to 
the first one in some norm. In other words, whether the solution to |(2.3)| exists and 
converges to the solution of the limit equation for as long as the General relativistic 
solution exists. 



3 The Postnewtonian limit as an initialization pro- 
cedure. 

The Theorem quoted above imply that solutions to Einstein's equations, (the system 
given by (2.20) to (2.22) ), can be splitted into slow and fast ones, and this is ruled 



by the choice of initial data. 

In the context of General Relativity the slow part of the solution contains infor- 
mation about the matter fields that are moving with velocity much lower than c, 
that is with its own time scale, while the fast part contains information essentially 
due to the arbitrariness of initial data and so is not related directly to the sources. 

In the context of symmetric-hyperbolic systems this means that on the Postnew- 
tonian approximations the highly oscillatory part (fast part) of the solution should 
be suppressed up to higher order of e. For example, if we initialize (choose initial 
data such that the time derivative of the solution is bounded independently of e 
at t = 0) up to order e the slow part appears at zeroth order and the fast one at 
first order. The degree of initialization needed depends on the problem at hand. In 
this Section we show that to get the quadrupole formula one necessary to initialize 
up to order e^. This guarantees that no contribution from fast solutions would be 
present in the gravitational radiation flux to that order. We then show that such 
initialization is possible. 

To know at which order of e we have to initialize, we estimate the contribution to 
the energy flux at X+ from both, sources and vacuum (pure initial data) solutions. 
Since we are only interested in the order of e at which these contributions appear. 



it is enough to estimate them in linearized gravity, for example see |]T2|. 
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To this end we keep the physical units in the Einstein equations (i.e k = G ^ 1), 
and we check the power at which e appears in the radiation due to the matter 
(quadrupole contribution, with zero initial data) and in the radiation due to the 
initial data (without sources). 

In this approximation 

and in the absence of matter the energy flux in the direction due to the initial 
data can be calculated as 

where the dot means the time derivative. In the presence of matter, with initial data 
^Qfiu = Sg^iu = 0, we can use the same formula as above because at large distance 
from the bodies the matter waves can be considered as plane waves. By |(2.19 
6gij = e'^hij, the flux of energy becomes 

1 



Thus choosing the initial data hij = e'^fij, the flux becomes ct°^ = 0{e'^°''^^), 
on the other side the flux due to the matter is calculated as ct°^ = 0{e^) (see the 
Appendix). We conclude that initializing up to order the fast part of the solution 
will contribute at a higher order of e to the flux. It means that the initial data will 
be choosen as 

u\t=o = uo + eul + + £^U3 + 0{e'^). 

This is equivalent to demand that the first four time derivatives of the solution 
stay bounded in t = when e —>■ (c.f. [§]). 

Besides an hyperbolic system Eintein's equations also include a constraint sys- 
tem, that is, equations which relate the otherwise arbitrary initial data fields between 
each other. Thus it must be chequed that these constraints are consistent with the 
initialization procedure. To do that we add to the equations arising from the initial- 
ization the constraint equations and show that in the resulting system, the hierarchy 
of elliptic systems admits solutions with the correct asymptotic fall off. 

Defining the pseudotensor /i"* 

(3.1) V^?'^* = v^e'^^ + ^'^"^ 

and from the definition of r"^c given by |(2.25)| , we obtain 



Choosing the gauge condition r'^^^ = and r°'''b = up to order and demanding 
the boundedness of the first four time derivatives we get that Pi""^ and /ij"*, i = 
0, 1, 2, 3 should satisfy at t = 0, the following differential equations: 



9 



Ao/io'^'' = 0, 

5ePo"' = 0, 

Ao/ii'^' = 0, 
Aoipi^'^ + d^^No'^) = 0, 
(3.2) Ao(p2"' + ^('^iVi'')) = 0, 

AoV* = 2Si'^^ - Fi"^^ + S^'^iVi^) - 2e"^Fi, 

where Aq = e"'^dadb. These equations are consistent with the constraint up to order 
6^, for they have been included in the set. 



Lemma 3.1 The equation system for the initial data given by \( 3. 2J\ admits a solu- 
tion. 

Proof: From the first five equations, we can choose 

(3.3) po"' = 0, V, = 0, ri"^ = 0, pi"" = -d^'^No'^ 
and 

(3.4) Pa"'' = -d^''Ni''l 

Therefore, in order to guarantee the existence of solutions, it is sufficient to prove 
that the sources in the three last equations in |(3.2)| have the decay 0{\) (c.f.[jl3[ 



Due to the choice made in PTT)! , the functions Fo"^ Fq, Fq^^, Fi, Fi" and Fi"'' 
are all linear functions of daiUoNo^) . Similarly Fq and Fq"'' depends linearly on 
dtdaiUoNf)') . But the vector fields Nq"" and Ni", and the function Uq satisfy the 
following differential equations: 

AoiVo'^ + 9"9eiVo^ = 4Jo^ 

Aof/o = Po- 

Thus, because of the compactness of the support of the sources, the standard 
theory of elliptic systems in R" (c.f. [0]) assures that Uq, Nq°- and Ni"-, and their 
time and spatial derivatives decay as C(^) and C(^) respectively. Hence we con- 
clude that the sources decay fast enough as to ensure existence of solutions for h2°'^, 
hs"^, and pg"* + ^("A/'s^). 

In order to have a well defined Ps"*, it remains to prove that N2"' exists. Note 
that the gauge condition implies that N2°' must satisfy 



10 



AoA^a'^ + d'^d.N^'' = 4J2" - SFi*^ - ANo^'dad'^Uo + 2Fo^\, 

and that the source in this equation has a fast enough decay. Thereby the existence 
of N2°' and so of ps"* is assured. □ 

This Lemma proves that we can initiahze at least up to order e^, hence by 
Theorem |2.1j we can assure that the fast solution appears in a higher order than 



Appendix: The evolution of the initiahzed data 

In order to exemplify the initialization procedure, we initialize up to order e and 
solve explicitly the equation for and v^, choosing the initial data uo and ui 

as in 1(3. 3)| . In fact this choice gives that u becomes a slow solution suppressing the 
fast one to a higher order of e. 



Claim: Given the initial data 



m(0, x) 



-di^^ Nb) 





(x) + 0{e' 



where N°'{x) is the solution of the elliptic equation ( 3. 6\ for t 
(2.21 \ and \( 2. 22\ can he written as 



0, the solution of 



with 



and 







0{e 



2\ 



pab ^ / -iK'^ ."''"^ dV + a^^iy") 

i|x-x'|<|- |X — X'l 

ff . I ' rf(t-.|x-x1.x')^3,, 

J|x-x'|<-^ X — X' 

I I — £ I I 

where Vf = -{Ad'-'^N^^ + 2e^S'''') and Tf = d'-'^N^^ + 2S''\ p and J" have compact 
support in R^. Furthermore, the asymptotic behavior of u is 0{-). 
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Proof: Since mq = and thereby K"'daUQ = 0, Theorem |2.1| assures that 
u{t, X, e) = u°{t, x) + e[v{t, x) + u{t, x, e)] + C(e^). 



By applying |(2.4)| and |(2.5)| on the equations system given by |(2.21)| and |(2.22) 
we obtain that and v will satisfy 



(3.5) K^daU^ = 0. 

At this order of e the constraint equations and the gauge become 



AU 

roc"' 



P 


4J^ 



(3.6) d.d^N^ + d^d^N'' - 

where h"-'' is a pseudotensor defined in |(3.1)| . Choosing 

/ _5(«iv^) \ 
' 



the system |(3.5)| becomes 
(3.7) 



u 



uo(x) 





0, 



which is trivially satisfied by u^{t,x) = 0. 

Considering now the equation for u, namely 



(3.8) 



with 



e 

m(0,x) = ui{x) -t;°(0,x) = 
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we obtain the following system 



c~ab 



(3.9) 

(3.10) r, + -d,p''' =0 

and the constraint and the gauge become 



p + -d'^r 

e 

r c + -OcP 



2Sab ^ Qiaj^b) 



(3.11) 



AU 

dcPT 



-2r 



ab 



4J^ 



Since p"'' = p'f — v""^ and Tc"^ = tic°*, instead of |(3.9)| we consider the differential 
equations 



5c„ ab 



(3.12) 



2gab 

e 



The equations system |(3.12)| is clearly consistent with the gauge |(3.11)| . Taking 
the time derivative to the system above, we obtain 



ab 



ab 



(3.13) 



2eTf 



and the solution of the second equation in |(3.13)| can be written as 



hf = 2e 



<'|<i 

I F 



vfit- 


- e 


X — 


x' 


,x') 




X 


-x' 





Because of the decay of the source, when |x| » |x'|, we write |x— x'| 
where ^ = Hence in this first approximation we have 



ab 



2e 

r J|x-x'|< 



/ rf (t-£r + en-x',x')dV, 

Jlx-x'Ki 



X — n-x 



similarly we calculate pf'^.D 

Following this procedure we can generate the solution that has as initial data 
1(3.2)1 and the functions v?{t, x) and u^{t, x) obey similar equations than vP{t, x) and 
u^{t, x) do. 
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Remark 3.1 In the calculus of the flux of energy due to the matter fields, we use 
the pseudo-tensor field calculated above instead of the spatial components of the 
perturbation of the metric. The difference between them can be calculated using 
equations \(2.19)\ and \(2.25)\ and because they have the same expansion in e they 
can be considered as equal. 
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